Quasi-2D Confinement of a BEC in a Combined Optical and Magnetic Potential 
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We have added an optical potential to a conventional Time-averaged Orbiting Potential (TOP) 
trap to create a highly anisotropic hybrid trap for ultracold atoms. Axial confinement is provided 
by the optical potential; the maximum frequency currently obtainable in this direction is 2.2 kHz 
for rubidium. The radial confinement is independently controlled by the magnetic trap and can be a 
factor of 700 times smaller than in the axial direction. This large anisotropy is more than sufficient 
to confine condensates with ~ 10 5 atoms in a Quasi-2D (Q2D) regime, and we have verified this 
by measuring a change in the free expansion of the condensate; our results agree with a variational 
model. 
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I. INTRODUCTION 

Two-dimensional systems are of great interest in 
condensed matter physics in general, and they have 
some special properties Q,0 A two-dimensional Bose 
gas confined in a uniform potential does not undergo 
Bose-Einstein condensation (BEC); instead there is a 
Berezinskii-Kosterlitz-Thouless (KT) transition: a topo- 
logical phase transition mediated by the spontaneous for- 
mation of vortex pairs which leads to a system that does 
not have long-range order but which is nevertheless su- 
perfluid. Many features of the KT transition have been 
observed in experiments with films of superfiuid He on 
surfaces 03; (see [j| for some recent results). In ex- 
periments with a thin layer of spin-polarised atomic hy- 
drogen gas on a liquid helium surface @, evidence was 
found for the existence of a quasi-condensate: a system 
that only possesses local phase coherence, in contrast to 
a pure condensate that has a global phase. 

The recent advent of laser cooling and subsequent 
Bose-condensation of the alkali metal atoms |7| has 
opened up a new avenue of investigation: while these 
bosonic gases are initially cooled to quantum degeneracy 
as three dimensional clouds, the addition of novel dipole 
force potentials can make these systems two dimensional 
to varying degrees. A condensate can be called quasi- 
two dimensional (Q2D) when the energy level spacing 
in one dimension exceeds the interaction energy between 
the atoms. In this case the particles obey 2D statistics 
but interact in the same way as in a three-dimensional 
system. 

The first Q2D condensate was made with sodium 
atoms that were condensed in a magnetic trap and then 
loaded into an attractive optical dipole trap The 
dipole trap was made by tightly focusing a laser beam 
along one direction with a cylindrical lens to achieve trap 
aspect ratios of 79. In this experiment the number of 
atoms, and hence the interaction energy, was reduced to 
enter the Q2D regime. A Q2D condensate of two thou- 
sand caesium atoms has also been made in a surface wave 
trap with an aspect ratio of 50 ■ In our experiment we 
use two sheets of blue-detuned light to increase the axial 



confinement of our magnetic trap, this approach allows 
anisotropies in excess of 700 to be achieved as the ra- 
dial confinement is provided independently by the mag- 
netic trap, which in turn allows more atoms to be loaded 
into the Q2D regime. The combination of our magnetic 
Time-averaged Orbiting Potential (TOP) trap and dipole 
potential creates a flexible system with the ability to ro- 
tate the condensate, so that vortex nucleation can also 
be studied. 

Q2D condensates have many interesting properties 
that are not present in 3D: the effect of phase fluctua- 
tions has been examined theoretically in |1C|: the possi- 
bility of vortex production as a mechanism for the de- 
cay of a quadrupole mode is examined in Indeed 
while it has been shown that in a harmonically trapped 
non-interacting gas Bose-Einstein condensation can oc- 
cur , there is still a question of when an interacting 
2D system supports a KT rather than a BEC transition. 

Our paper is structured as follows: in the second sec- 
tion we describe the design of our optical potential, in 
the third section we discuss its realisation and some ini- 
tial calibration measurements, and in the final section we 
present our measurements of the condensate in free ex- 
pansion which confirm it has been trapped in the Q2D 
regime. 



II. DIPOLE TRAP DESIGN 

A. Optical Potential 

In order to produce strong axial confinement, we trap 
the atoms in the nodal plane at the focus of a Hermite- 
Gaussian TEMoi-like mode. To produce this light field 
the output from a laser at wavelength A = 532 nm (Co- 
herent Verdi) is passed through a phase-plate, which im- 
prints the top half of the beam with a tt phase shift rel- 
ative to the bottom half. The phase-plate was manufac- 
tured in house, and coated in such a way as to balance 
the reflective losses from the upper and lower surfaces. 
After focusing with a cylindrical lens with focal length 
/ = 160 mm in the z direction, the intensity pattern at 
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the focus resembles a TEMoi mode to a first approxima- 
tion. 

We derive the exact expression for the intensity distri- 
bution below and also discuss any shortcomings together 
with the effects of various types of misalignment. The 
axis convention used in this paper is explained in Fig. ^ 
and the coordinate axes are denoted by capital letters for 
the input plane (at the position of the cylindrical lens) 
whereas they are denoted by small letters in the focal 
plane, where the atoms are trapped. The intensity dis- 
tribution of the initial beam is given by 
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where a y ,a z are the beam waists in Y and Z direction, 
respectively, and P is the power in the laser beam. The 
amplitude function (f>(y, z) in the focal plane is given by 
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where 4>{y, Z) = cp{Y, Z) = y/I init (Y,Z) j(Z - Z ) is 
the amplitude function in the input plane. The function 
j(Z - Z Q ) = 1 for Z < Z and j(Z - Z Q ) = -1 for 
Z > Zq. It represents the action of the phase-plate which, 
as we shall assume in the following, is aligned with the 
center of the Gaussian beam for Z = Zq- The cylindrical 
lens performs a 1-d Fourier transform and k — kz/f, 
k = 2tt/X. After some calculations we obtain for the 
intensity in the focal plane 
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where the imaginary error-function erfi(z) is defined in 
terms of the conventional error function erf(z) evaluated 
for a purely imaginary argument: erfi(z) = — i crf(iz). 
To first order erfi(z) as (2/y / 7r)z + 0(z 3 ). Inserting this 
approximation to erfi. into the above equation and ex- 
panding it up to second order in z we obtain 



Ifoc(z) = 
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Provided that A ^> Aps UM' where A is the detuning of 
the dipole trap from resonance with the atomic transition 
and Aps is the fine structure splitting of the transition, 
the potential energy of an atom in the dipole trap is 
related to the intensity of the light field by 



E^ p (r) = ^I/frl 
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where c is the speed of light in vacuo, T and are the 
linewidth and frequency of the atomic transition, and 
I(r) is the intensity of the dipole trapping beam. In our 



experiment the atomic transition is the 2 <Si/2 -^> 2 P\ 12,3/2 
in atomic Rb 87 which has a wavelength of 780 nm. From 
Eas. l)4l5| l we obtain for the axial trap frequency Vt, for 
atoms of mass m, 
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Going back to Eq.® we shall define the beam waists in 
the focal plane as w y — a y and w z = 2f/(ka z ). The 
latter relation is exactly that between the input and the 
output waists of a Gaussian beam. A Gaussian beam of 
any order transforms into different sized Gaussian beam 
of the same order through the action of a lens. In our 
case the phase jump introduced by the phase-plate causes 
a zeroth order Gaussian to transform into something 
closely resembling a first order Gaussian beam (TEM01) 
rather than a zeroth order. Using the expressions for the 
beam waists in the focal plane the intensity profile © 
can be written as 
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The axial trap frequency £1 can then be expressed in 
terms of these waists as 



n 2 = 
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In order to determine the trap frequency f2 for our exper- 
imental parameters we can measure the TEMoo profile in 
the input plane and then use Eq.© to calculate the re- 
sulting trap frequencies of the optical potential in the 
focal plane. Alternatively, we can measure the distance 
d z between the two intensity peaks in the focal plane 
and use the relation d z = 1.8483 w z to determine fl from 
Eq.(jHJ. The beam intensity profile, given by Eq.© dif- 
fers somewhat from that of a perfect TEM01 which is 
given by 
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Substituting erfi(z) w 2/(y / 7r)z into Eq.Q) we obtain 
Eq.© with the initial multiplying constant smaller by 
a factor of w. This approximation is only valid for small 
z but it indicates that in the central region, where we 
want to confine the atoms, the intensity of the focused 
beam is reduced by a factor of tt with respect to that of 
a pure TEM01 beam of the same power P. The reason 
is that much power is lost in the very broad wings of the 
focused beam which are much larger than those of a pure 
TEM01. This can be seen in Fig|2Ji, where the solid line 
shows the actual potential given by Eq.0 and the dot- 
ted line shows the potential of a TEM 01 given by Eq.@ 
for comparison. 



3 



The initial amplitude distribution has a sudden jump 
in the middle where the amplitude goes from negative to 
positive. The focusing lens performs a Fourier transform 
and in order to resolve the sudden change in amplitude 
many higher spatial frequencies, represented by the wings 
of the intensity profile in the focus, are needed. 

Combining the magnetic with the optical potential we 
find that the central position of the combined potential, 
in the axial direction, is given by 

*c(t) = ^5-, (10) 

where z c is the position of the combined potential rela- 
tive to the center of the magnetic trap, zi is the position 
of the center of the optical potential, and loi and u rn are 
the angular frequencies of the optical and magnetic po- 
tentials respectively. 

B. Beam Shaping and Misalignment Errors 

A schematic diagram of the experimental setup is 
shown in Fig. 01 T° achieve a high axial confinement fre- 
quency, a tight focus is required in the z direction, as such 
the beam is expanded in this direction prior to focusing. 
A suitable choice of beam width is also required in the y 
direction to determine the area in the xy plane through- 
out which Q2D confinement can be realized. The y and 
z beam widths are not the same, and so some method 
of asymmetric beam expansion is required: for this ex- 
periment two orthogonal, cylindrical telescopes are used. 
To achieve a beam focus of 7.2 fj,m with a final lens of 
focal length f=160 mm, the beam waist in the z direction 
before focusing must be 3.8 mm. The beam waist in the 
y direction is chosen to be 410 /xm. From Eq. © the 
resulting axial trap frequency is w 2.7kHz. 

The potential height and frequencies are decreased by 
misaligning the phase-plate with respect to the incoming 
Gaussian beam. Two types of misalignment can occur. 
The first one is displacement of the phase step on the 
phase-plate from the center of the Gaussian beam. In 
this case Z , defined below Eq.(0|, is not equal to zero. 
As a result the intensity dip in the central region is fills up 
for large values of the displacement parameter £ = Z j a z . 
Fig- Eb shows a series of potential plots for various values 
of £. The central dip has completely disappeared for £ = 
1.5. The resulting change in trap frequencies is shown in 
Fig. 01 (solid line) which plots normalised by its value 
when £ = 0. The dotted line plots the potential height 
normalised by its value for £ = 0: A pot = (I m ax(£) ~ 

The effect of rotational misalignments between the in- 
put beam, the phase-plate and the final cylindrical lens 
have been calculated numerically for our experimental 
parameters (see Fig.[SJ). The results show that the optical 
potential at the beam focus is relatively insensitive to the 
precise rotational alignment of the phase-plate; however, 



the angle that the elongated input beam makes with the 
cylindrical lens is much more critical, as misalignment de- 
creases the axial confinement and introduces some radial 
confinement. In our experiment we aim to minimise the 
radial confinement and maximise the axial confinement 
provided by the optical potential. 



III. EXPERIMENTAL REALISATION 
A. Alignment 

Sub-micron actuators give independent control over 
the focus of the optical potential in three dimensions. 
The x and z position of the focus can be controlled by 
moving the final cylindrical lens, and the y position is 
changed by moving the position of a mirror which de- 
flects the beam through 90 degrees (see Fig. [3Jl. Inten- 
sity control is provided by an Electro-Optic Modulator 
(EOM). 

The size and depth of the optical potential are too 
small to provide any noticeable change to the Magneto- 
Optical Trap (MOT) used in the first stage of our ex- 
periment, and so initial alignment was carried out by 
making the beam from a 397 nm semiconductor diode 
laser co-linear with the dipole trapping beam. Photons 
at this frequency are sufficiently energetic to ionize Rb 
atoms and deplete the MOT, providing a clear signature 
for alignment. Subsequent coarse alignment was carried 
out by lining up the MOT with the entry and exit reflec- 
tions of the trapping beam on the glass cell. This allows 
the beam to be aligned to within 200 /im, in x and y, of 
the magnetic trap center. 

For the next stage of alignment the atoms were sup- 
ported against gravity on top of the laser beam. The 
power in the trapping beam was then reduced, until the 
atoms were only just supported, and the beam moved 
in the xy plane: as the beam focus is moved closer to 
the atoms, less power is required to stop the atoms from 
falling. This is an iterative process, which is accurate to 
~ 30 fim in x and y. 

The most critical direction for alignment in our exper- 
iment is the z direction, which needs to be aligned to 
within a few microns to load atoms from the magnetic 
trapping potential into the dipole trap. Alignment in 
this direction exploits a delay that is required between 
the switch off of the magnetic and optical potentials. 
The dipole trap can be switched off by the EOM in a 
few microseconds, however the magnetic trap requires 
almost 500 /is before the current in the quadrupole coils 
decays. This can create a problem when looking at the 
expansion of condensates that are tightly confined by the 
optical potential: simulations show that when suddenly 
released, the expansion of the condensate can be slowed 
or even temporarily reversed by the residual magnetic 
fields. The solution is to introduce a delay of 500 /is af- 
ter the magnetic field begins to decay, before switching 
off the optical potential. The equilibrium position of the 
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atoms in the combined potential is not quite the same as 
the equilibrium position in the purely optical potential: 
the result is that the atoms move in the optical poten- 
tial during the delay. The velocity they acquire during 
this period is sufficient after 15 ms of free expansion to 
allow the start position to be calculated, and is accurate 
to ~ 1 /im. 

In Fig.[5]the position of the dipole trap beam is scanned 
in the z direction, while the magnetic trap, which de- 
fines the initial position of the atoms, remains fixed. The 
dipole trap intensity is then ramped on and the position 
of the atoms observed after free expansion. Fig. [fj{c) 
shows the pattern observed and expected for a TEM i 
potential, as shown in Fig.[|Jd). The magnetic and dipole 
traps are only aligned in the central region of the plot in 
Fig. ini^c), this coincides with a sharp increase in the ax- 
ial expansion of the condensate. The beam waist can be 
accurately measured from the width of the central re- 
gion. This method can also be used to profile the beam 
at various points along its optical axis, and to partially 
reconstruct the intensity profile. Fig. Ela) shows a scan 
further away from the beam focus, and Fig. IHIb) is an 
intensity pattern that produces the observed deflection. 



B. Loading the Trap 

Our experiment uses evaporative cooling in a TOP trap 
to create a Bose-condensate that contains ~ 10 5 atoms 
of the Rb 87 isotope in the \F = 1, mp = —1) state. The 
initial oscillation frequencies of atoms in the magnetic 
trap are 61 Hz radially and 172.5 Hz along the z (ax- 
ial) direction. After a condensate has been formed it is 
adiabatically loaded into the combined potential: the in- 
tensity of the dipole force laser beam is ramped from zero 
to its maximum value in 100-300 ms (depending on the 
strength of the optical potential). This ensures that no 
dipole modes, or internal modes of the condensate, are 
excited during the loading phase. This was confirmed by 
an exact numerical simulation. The combined central po- 
sition changes most rapidly when the frequencies of the 
potentials are comparable (see Eg. ITU)) , at the beginning 
of the ramp. For the highest frequency dipole traps this 
means that the ramp needs to be extended to give the 
condensate a smooth transition and avoid unnecessary 
heating. 



C. Calibration 

Dipole frequency measurements in the radial and ax- 
ial directions are used to characterise the trap. For the 
highest frequency traps, the measurement is complicated 
because the TOP frequency (7 kHz) is comparable to the 
dipole frequency in the axial direction. Efforts to ob- 
serve the dipole oscillation with the TOP field on were 
unsuccessful for this reason. It is therefore necessary to 
turn off the TOP field and move the quadrupole center 



above the optical potential, to avoid atom loss through 
Majorana spin flips. In this geometry it is necessary to 
suddenly change the quadrupole gradient to excite the 
dipole mode. We have measured an axial frequency of 
2.2 kHz, as shown in Fig. [3 for a beam waist of 7.2 /im, 
determined from the deflection of atoms as previously 
explained. This agrees roughly with the expected value 
of 2.7 kHz for an ideal phase-plate profile, calculated us- 
ing Eq. when the rotational misalignment between the 
phase-plate and the final cylindrical lens is taken into 
account. A rotation by only 4 ± 1 degrees effectively 
changes the axial frequency to 1.94 =F 0.27 kHz. From the 
frequency measurement it is also possible to estimate the 
axial beam drift to be approximately 18 /im/hr. The ra- 
dial trap frequencies have also been measured by exciting 
dipole modes, and by measuring the aspect ratio of the 
condensate (defined as l y /l x -, where / is the width of the 
condensate extracted from a Gaussian fit) in the xy plane 
after free expansion. At high values of radial confinement 
from the magnetic trap, the dipole frequencies were as ex- 
pected for a purely magnetic potential; however, as the 
magnetic confinement was reduced it was observed that 
the cloud elongated along the beam direction (see Fig.[SJ), 
and when the dipole mode was excited it oscillated only 
along the propagation direction of the trapping beam. In 
this direction the dipole oscillation frequency agrees with 
the expected value for a purely magnetic potential, but 
across the beam there is some extra confinement provided 
by the optical potential. The optical frequency along the 
y direction was used as a fitting parameter for the data 
in Fig. |SJ for a fixed optical frequency a curve was cal- 
culated for all the various magnetic frequencies and the 
resulting curve compared to the data points. This pro- 
cedure was repeated to find the least squares minimum 
value for the optical potential. This was found to be 
equivalent to 26 Hz for an axial trapping frequency of 
1990 Hz, and is caused by the rotational misalignment 
described above. Indeed subsequently we have measured 
the tilt of the input beam relative to the final cylindrical 
lens and found a value of 4 ± 1 degrees, from Fig. [S] it 
can be seen that this theoretically gives between 19 and 
23 Hz of confinement. 

We also measured the axial size of the condensate as 
a function of trapping beam power. The condensate ex- 
panded for 15 ms before it was imaged. The results are 
shown in Fig. [5] and are in good agreement with the hy- 
drodynamic theory and the assumption that f2 cx V^P, as 
would be expected for a harmonic dipole trap (see Eq.[f)J). 
These measurements were made at a radial frequency of 
61 Hz and so the hydrodynamic expansion theory applies. 



IV. THE QUASI-2D REGIME 

A. Criterion for Q2D 

The condition fi < huj zi where /i is the chemical po- 
tential, can be used as a criterion for when a condensate 
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is in the Q2D regime. In this case the ground state en- 
ergy is smaller than the harmonic oscillator spacing and 
thus comparable to the harmonic oscillator ground state. 
When ~ 10 5 atoms have been loaded into our combined 
optical and magnetic potential the Q2D regime can be 
reached by either decreasing the number of atoms, or by 
adiabatically decreasing the radial confinement provided 
by the magnetic trap. In the following experiments we 
use the second approach, to provide more control. Along 
the z direction the condensate shape becomes very simi- 
lar to the Gaussian profile of an ideal gas. However, along 
the weakly confined x and y axes the condensate is char- 
acterised by the hydrodynamic parabolic shape. For the 
best description in terms of simple analytic functions it 
is therefore best to describe the condensate wavefunction 
as a hybrid of a parabolic and a Gaussian distrubution. 

B. Expansion of a Q2D gas 

In our experiments we observed the condensate ex- 
pansion in various regimes and the smooth transition 
from the hydrodynamic expansion characteristics to 
those of a quasi-2D gas, which essentially expands like 
an ideal non-interacting gas in the direction of tight con- 
finement. The results are shown in Fig. llOf a). where the 
open circles and filled circles are the data for traps with 
u} z /2ir = 1990 Hz and 960 Hz respectively. The expan- 
sion time is constant at 15 ms and the radial trap fre- 
quency is varied to explore the transition to Q2D. These 
curves demonstrate clearly the ideal-gas like behaviour 
for the axial expansion in this limit. We compared our 
findings to those of theoretical predictions which were de- 
rived from variational models [15j and found good agree- 
ment to our data. The prediction of the hybrid varia- 
tional model 0| , which is indistinguishable from that of 
the Gaussian variational model, is given by the solid line. 
The horizontal dashed lines indicate the expansion of the 
ideal gas, given by R z — ty/fiu z /m, where t is the ex- 
pansion time, and the dotted lines the expansion of the 
hydrodynamic gas. The 'ideal gas' and hydrodynamic 
models yield straight lines on the logarithmic plot. In 
contrast, the curve for the hybrid model follows the hy- 
drodynamic asymptote down to u) x ~ 20 Hz, correspond- 
ing to [i « TiLOz, where it bends and follows the 'ideal gas' 



line towards zero radial frequency. The two noticeably 
do not coincide at ui x w because of the residual opti- 
cal anisotropy in the radial plane (as we discussed before 
in the text describing Fig. 0, which was taken into ac- 
count in our variational model and is in good agreement 
with the experimental data. This transition from the hy- 
drodynamic to the 'ideal gas' asymptote gives conclusive 
evidence of the gas entering the Q2D regime. The release 
energy of the condensate, calculated from the size after 
expansion, is also shown in Fig. llOf b): this tends towards 
the ideal gas limit of Tilj z /A as uj x is reduced. This is the 
vertical kinetic zero-point energy; all that is available if 
the confining potential is suddenly switched off. 

V. CONCLUSIONS AND OUTLOOK 

In these experiments we have studied the properties 
of quasi-2D condensates. To obtain the required ge- 
ometry we increased the stiffness of the axial confine- 
ment by superimposing an optical potential on the exist- 
ing magnetic trapping potential. The resulting trap was 
characterised in detail and we obtained confinement fre- 
quencies close to the theoretically expected values. We 
then studied the properties of Q2D condensates including 
their chemical potential and expansion characteristics. 
Our experimental observations agreed with theoretical 
predictions to confirm that we successfully entered the 
Q2D regime. This was the first experiment to explore 
the hydrodynamic-Q2D transition by gradually increas- 
ing the trap anisotropy instead of throwing atoms away 
as in previous work |8|. In future experiments we hope 
to study the properties of vortices and vortex arrays in 
the Q2D regime and investigate the possibilities of ob- 
serving a Berezinskii-Kosterlitz-Thouless transition and 
studying its properties. 
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FIG. 1: TEMoi-like intensity distribution at the focus of the 
final cylindrical lens. Axis convention is shown: z is in the 
vertical direction, y is horizontal across the beam and x is 
along the direction of beam propagation, out of the page. 
The y and z directions are not to scale. The radial direction 
is to be interpreted as in the xy plane. 
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FIG. 2: (a) The intensity distribution in the focal plane for 
the actual experimental setup (solid line) and for the ideal 
TEMoi profile (dotted line) of the same total power, (b) The 
actual potential is plotted for various values of the displace- 
ment parameter: £ = 0.3 (solid line), £ = 1.0 (dashed line), 
£ = 1.5 (dotted line). 
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FIG. 3: Schematic drawing of optical layout showing the 
cylindrical telescopes, the n phase-plate and the final cylin- 
drical lens. Shaded arrows indicate the location of sub-micron 
actuators for positional control over the focus. 




FIG. 4: The oscillation frequency for rubidium atoms in the 
potential normalised by its value for £ = {£l/Q. max ) is plot- 
ted against the phase-plate displacement parameter £ (solid 
line). The potential height A po t decreases with increasing £ 
(dotted line). 
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FIG. 5: The effect of angular misalignment between the input 
beam and the final lens, calculated for w z —7.2 fim and P = 
4.5 W. (a) Plot shows the confinement in the y direction as a 
function of angular misalignment, (b) Plot shows decrease in 
axial trapping frequency as a function of misalignment, (c) 
Intensity distributions at the focus of the cylindrical lens for 
2, 5 and 10 degrees of input beam misalignment (left to right; 
in all the plots the y and z axes are not the same scale). 
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FIG. 6: The position of the dipole trap beam is scanned in the 
z direction, while the magnetic trap defining the initial posi- 
tion of the atoms, remains fixed. The dipole trap intensity is 
then ramped on and the position of the atoms observed after 
15 ms of free expansion. (a)/(c) show a scan far from/closer 
to the focus of the trapping beam, 'white squares' are exper- 
imental points and the 'solid line' is the deflection expected 
for the potential in (b)/(d). 
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FIG. 7: Axial frequency measurement with TOP field turned 
off. 'Black squares' are experimental points, 'solid line' is a 
sinusoidal fit to the data used to extract the frequency. 




FIG. 8: Aspect ratio of the condensate, defined as a y /a x , after 
20 ms free expansion. 'White triangles' and 'black squares' are 
the measured aspect ratios of a cloud released from a hybrid 
trap with / z =960 Hz and f z = 1990 Hz respectively. 'Solid 
line' and 'dashed line' are fits to the data using the varia- 
tional method; this allows the frequency of the optical poten- 
tial along y to be determined. For the data given / H = 12.0 Hz 
and f y =26.0 Hz just from the optical potential, as expected 
there is more confinement in the trap with the higher axial 
frequency. 
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FIG. 9: Axial condensate size after 15 ms of free expansion, 
as a function of power. The axial/radial trap frequency is 
960 Hz/61 Hz. 'Black squares' are the experimental data, 
'black line' is the hydrodynamic prediction provided that f z oc 
VP, as is expected for a harmonic potential. 
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FIG. 10: (a) The axial expansion of the condensate after 
15 ms of free expansion for various trap geometries and atom 
numbers. Solid lines indicate theoretical variational predic- 
tions, dashed lines indicate the ideal gas limit and dotted lines 
the hydrodynamic limit. The data are taken for traps with 
u) z /2tt = 1990 Hz (open circles) and 960 Hz (filled circles). 
The atom numbers are 8 x 10 4 and 1.1 x 10 5 , respectively, 
(b) The release energy of the condensate derived from the ex- 
pansion measurements in (a). The energy tends towards the 
vertical zero-point kinetic energy as lo x is reduced. 



